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Abstract: We study the axial anomaly denned on a finite-size lattice by using a 
Dirac operator which obeys the Ginsparg- Wilson relation. When the gauge group 
Ph. is U(l), we show that the basic structure of axial anomaly on the infinite lattice, 

which can be deduced by a cohomological analysis, persists even on (sufficiently 
>► ! large) finite-size lattices. For non-abelian gauge groups, we propose a conjecture on 

a possible form of axial anomaly on the infinite lattice, which holds to all orders 
in perturbation theory. With this conjecture, we show that a structure of the axial 
anomaly on finite-size lattices is again basically identical to that on the infinite lattice. 
Our analysis with the Ginsparg- Wilson Dirac operator indicates that, in appropriate 
frameworks, the basic structure of axial anomaly is quite robust and it persists even 
in a system with finite ultraviolet and infrared cutoffs. 
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1. Introduction 

In ref . [1] , Liischer pointed out that a cohomological analysis can be used to determine 
a basic structure of the axial anomaly in abelian gauge theories with finite lattice 
spacings. This work paved a way to study the axial anomaly in a system with a finite 
ultraviolet cutoff and then the technique was applied for various cases [2, 3, 4, 5, 6]. 
The crucial properties which make this analysis possible are the locality, the gauge 
invariance and a topological property of the axial anomaly. The axial anomaly de- 
fined by the gauge covariant Dirac operator [7, 8] which satisfies the Ginsparg- Wilson 
relation [9], especially the overlap-Dirac operator [8], in fact possesses the required 
properties [10, 11, 12, 13, 14]. A further elaborate analysis with this recognition 
finally led to a non-perturbative construction of anomaly-free abelian chiral gauge 
theories on the lattice [15]. 

The cohomological analysis, however, is limited to the case of a lattice with an 
infinite size. A direct cohomological analysis for finite-size lattices is not feasible 
because (i) The analysis is based on the lattice Poincare lemma [1], which is a lattice 
analogue of the Poincare lemma being valid for R d . When the topology of the lattice 
is non-trivial (as is the case for the periodic lattice), one expects a non-trivial d- 
cohomology on the lattice, (ii) The cohomology relevant to an analysis of axial 
anomaly is a local cohomology, for which the concept of the locality is vital. The 
meaning of the locality, however, is not clear on a lattice with a finite size because a 
Dirac operator which obeys the Ginsparg- Wilson relation has to have exponentially 
decaying tails [16, 17]. 

In this paper, we study the axial anomaly defined on a finite-size lattice by 
using the Ginsparg- Wilson Dirac operator. This analysis provides an approach to 
the axial anomaly in a system with ultraviolet and infrared cutoffs. As already 
noted, a direct generalization of the technique of ref. [1] is not feasible. Instead, 
we point out that it is possible to determine the structure of axial anomaly using 
an argument similar to that of ref. [15] at least in abelian gauge theories. For non- 
abelian theories, we propose a conjecture on a possible form of axial anomaly on the 
infinite lattice, which is correct within perturbation theory. Under this conjecture, 
a similar argument can be applied to non-abelian cases too. These results indicate 
that the structure of axial anomaly is quite robust even with ultraviolet and infrared 
cutoffs in appropriate formulations (in the present formulation based on the 

Ginsparg- Wilson relation). We consider an even-dimensional lattice r whose size 
is L, r = {x G Z d | < x^ < L}, and the gauge field U(x,fi) G G (G is the 
gauge group) is assumed to be periodic on r, U(x + Lu, /i) = U(x, /i). 1 The lattice 
spacing a is set to be unity, except when the classical continuum limit is considered. 



ji denotes the unit vector in direction \i. 
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2. Preliminaries 



The axial anomaly for the Ginsparg- Wilson Dirac operator is defined by (see for 
example refs. [11, 12] for the background) 



A(x) = tr7 d+ i 



1 - ^D(x,x) 



(2.1) 



The kernel of the Dirac operator D(x,y) satisfies the Ginsparg- Wilson relation 

ld+1 D(x, y) + D{x, y) ld+1 = D ld+1 D(x, y). (2.2) 
The salient feature of A(x) is a lattice analogue of the analytic index theorem [10] 

J2A(x)=n + -n., (2.3) 

which follows from the algebraic relation (2.2) alone; here n + (n_) is the number of 
zero-modes of 7d+i-D with the positive (negative) chirality. The index theorem (2.3) 
implies that the Dirac operator cannot be a smooth function of the gauge field in 
general, because the configuration space of lattice gauge field is arcwise connected 
and, barring a possibility that n + — is constant for all configurations, the inte- 
ger n + — n_ jumps at certain points in the configuration space. A sufficient condition 
for the smoothness of the overlap-Dirac operator [8] is the admissibility [13, 14] 

|| 1 — U{x, v)\\ < e, for all x, v, (2.4) 

where U (x, fi, v) is the plaquette variable and e is a constant smaller than (2 — 
\f2) / d(d - I) [IA] 2 After imposing this admissibility, the space of allowed gauge 
field configurations may have non-trivial topology. This condition also guarantees 
the locality of the operator [13, 14] 

\\D(x,y)\\ < C(l + \\x - y\\ p )e-^-^ /e , (2.5) 

where C and p are constants and g is a localization range of the Dirac operator. In 
addition to the gauge covariance and the locality of the Dirac operator, we assume 
that it has the same transformation law as the standard Wilson-Dirac operator under 
discrete symmetries of the lattice (rotations, reflections, etc.). In particular, we 
require the translational invariance, i.e., D(x, y) is identical to D(x + z, y + z) if the 
gauge field is shifted at the same time U(x,fi) — > U(x + z, [/,). 

Suppose that we have constructed a Dirac operator on a lattice with the size L. 
When L — > oo, D(x,y) is promoted to a Dirac operator on the infinite lattice 
D(x,y) — * D°°(x,y). This operator also obeys the Ginsparg- Wilson relation 

ld+l D°°{x, y) + D°°{x, y) ld+l = D°° ld+1 D°°(x, y). (2.6) 



2 When the mass parameter m in the overlap-Dirac operator is unity, \m\ = 1. 
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In what follows, when we compare objects on the finite lattice r and on the infinite 
lattice, we always take repeated copies of a configuration of the gauge field on r as 
the gauge-field configuration on the infinite lattice. D(x,x) on r and D°°(x,x) with 
the argument x restricted to r may somewhat differ because they have exponentially 
decaying tails. However we assume that this "finite size correction" is exponentially 
small, 

\\D(x, x) - D°°(x, x) || < KL u e~ L/e , for ief, (2.7) 

where k and v are constants. The overlap-Dirac operator [8] possesses all the required 
properties we assumed above. 3 

Now, on the infinite lattice, the axial anomaly is given by 

1 



A°°(x) = tr 7d+1 
This is a topological field in the sense that 



2 -£>°°(x,x) 



(2.8) 



£ 5A°°(x) = 0, (2.9) 

where S denotes a local variation of the gauge field. This property can be shown 
from the Ginsparg- Wilson relation (2.6) (see ref. [2] for example). A°°(x) is thus a 
local topological gauge invariant pseudoscalar field. 4 When the gauge group is U(l), 
we can then apply the cohomo logical analysis [1, 2] to this quantity. The result is 5 

A 0O (x) = q(x) + dlk?(x), (2.11) 

where k^(x) is a local gauge invariant axial vector current (which is translational 
invariant) and the topological density q(x) is given by 

(470^/2)! 

xF M^d/i( x + fa + ^ H h frd/2-1 + "d/2-i), (2.12) 

with an integer Af. The abelian field strength is defined by 6 

F lu/ (x)=hiU{x, f ji,u), -7T < \ F^{x) < ir. (2.13) 



3 For the overlap-Dirac operator, whose basic building block is the Wilson-Dirac operator, one can 
show the relation D(x, y) — J2 ne zd D°° (x , y + Lri) . We thank Yoshio Kikukawa and Martin Luscher 
for clarifying this point. From this relation and the locality (2.5), one obtains the bound (2.7). 

4 A field 4>{x) is termed local, when its dependence on the gauge field at a point y is exponentially 
suppressed as \\x — y\\ — > co. For a more precise definition, sec ref. [1]. 
9 M and d* 



5 d ll and d* denote the forward and the backward difference operators respectively: 



d li f(x)=f(x + p)-f(x), d*f(x) = f(x) - f(x - p,). (2.10) 

6 For the cohomological argument to apply the constant e in eq. (2.4) has to be smaller than f 
and \F^(x)/i\ < tt/3 [1]. 
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Strictly speaking, the cohomological analysis alone admits a more general form of q(x) 
than eq. (2.12); for example, P^F^x) with anti-symmetric constants (5^ is also 
possible. However, since A°°(x) is a pseudoscalar under lattice rotations and re- 
flections, one infers that it must be proportional to the Levi-Civita symbol. Also 
the numerical coefficient in eq. (2.12) is left undermined in the cohomological anal- 
ysis. We can however use a matching with the result in the classical continuum 
limit; the integer Af is given by a sum of chiral charges of massless degrees of free- 
dom [18, 19, 20, 21, 22, 23, 24]. 

Note that eq. (2.11) is a statement for finite lattice spacings. Eq. (2.11) states 
that, even when the lattice spacing is finite, the main part of the axial anomaly 
is given by the topological density q(x) which has a quite analogous form to the 
continuum counterpart. On the other hand, the total divergence term d*k°f(x) 
represents "lattice artifacts" in the axial anomaly which depend on the details of 
the Dirac operator adopted. Our aim in this paper is to show or argue that the 
structure represented by eqs. (2.11) and (2.12) persists even on finite-size lattices 
and for general gauge groups G. 

3. Abelian gauge theory G = U(l) 

For the axial anomaly defined on a finite lattice (2.1), a direct cohomological analysis 
is not feasible. Nevertheless, we can show the following 

Theorem 3.1 When G = U(l) ; if the lattice is sufficiently large compared to the 
localization range g of the Dirac operator, say L/g> n, 

A(x)=q(x)+d;k,(x), (3.1) 

where k^(x) is a gauge invariant periodic current on r. The current k^(x) moreover 
satisfies the bound 

IM*)-*~(*)l<Ki^e~ L/e » (3-2) 

with constants K\ and v\. 

We emphasize that, for a sufficiently large L, eq. (3.1) is an exact statement for 
the axial anomaly A(x). Eq. (3.2) shows that the current k^(x) differs from the 
local current k^(x) defined on the infinite lattice only by an exponentially small 
amount. Hence, when the lattice size becomes large compared to g and thus when 
the concept of the locality becomes meaningful, the current k^(x) can be regarded 
as a local current. In this way, eq. (3.1) shows that the structure of axial anomaly 
on finite-size lattices is basically identical to that on the infinite lattice (2.11). The 
validity of this theorem has been argued intuitively by Chiu [25]. 
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Proof. The configuration space of the gauge fields allowed by the admissibility (2.4) 
consists of many components. Each component is uniquely characterized [15] by the 
magnetic flux 



1 L - x 



"V = 7T~ J2 F liv (x + sfi + tO), (3.3) 

zm s,t=0 



which is an integer. For a configuration with the magnetic flux m^, from eq. (2.12) 
one has [26] 

M(-l) d ' 2 



E A (x) — ^ q(x) — 2d/2^^2)! e ^ lUl '"M2 U d/2 m f J -l^i m ^2 ' ' ' m Md/2^d/2 

= an integer, (3.4) 



where the first equality follows from the translational invariance of k°f(x) (namely, 
k^(x) is a periodic current on r, when the gauge field is periodic). 7 Combined with 
the index theorem (2.3), we see that J2xer A(x) — J2xer A°°(x) is an integer. This 
integer is however bounded by an exponentially small quantity: From the assumed 
property (2.7), one infers that 

£ [A(x) - A°°(x)\ < K 2 L U2 e- L / s . (3.5) 

Therefore, when L is greater than some multiple of g, one has 

Y,[A(x)-A°°(x)}=0. (3.6) 

xer 

For this, we can apply the following 

Lemma 3.1 For a periodic field c(x) on r satisfying 

E c ( x ) = 0. ( 3 - 7 ) 
xer 

there exists a periodic current b^(x) which is given by a sum of c(y), the precise 
meaning of which is given in eq. (3.9) below, such that 

d*bJx) = c(x), \bJx)\ < 2Lmax|c(a;)|. (3.8) 
" xer 

Applying this lemma to eq. (3.6), we see that there exists a gauge invariant periodic 
current Ak^(x) such that A(x) — A°°(x) = d^Ak^x). This field is exponentially 
small, \Ak^(x)\ < « 1 L 1/1 e _L/ ' e , thus k^(x) = k^(x) + Ak^(x) which proves the theo- 
rem. □ 



7 Eq. (2.3) and theorem (3.1) show that the index is given by the combination (3.4) in terms of 
the magnetic flux. For the overlap-Dirac operator, this relation has been verified numerically for 
d = 2 and d = 4 [27, 28] and proven analytically for d — 2 [28]. 
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The assertions of the lemma immediately follow from the explicit construction 
of (though this is not unique): 

^ L-l L-l 

h ^ x ) = Jd^lYl ■•• S c(xi,...,x M _i,^, ...,y d ) 

x +1 L_1 L ~ l 

S ••• S c(xi,...,x M _i,j/ M ,...,j/ d ). (3.9) 



Note that since b^(x) is given by a sum of the field c(x), b^x) is gauge invariant if 
so is c(x). 



4. Non-abelian cases 

For general gauge groups G, a cohomological argument in a non-perturbative level is 
not known even on the infinite lattice. Thus we propose a conjecture on a possible 
form of A°°(x): 

Conjecture 4.1 For general G, 

A°°(x) = q(x) + d;k~(x), (4.1) 

where k?f(x) is a local gauge invariant axial vector current (which is translational in- 
variant) and the topological density q(x) is given by Liischer's topological density [29] 
and its higher dimensional extensions. 

The explicit expression of Liischer's topological density is known only for d = 2 and 
for d = 4: In our context, it is given by Af times eq. (32) of ref. [29]. We simply 
assume that the construction can be pursued for higher dimensional cases. 8 The 
construction of ref. [29] does not provide a pseudoscalar q(x). However, we may 
always enforce this pseudoscalar property by taking average over lattice symmetries; 
we assume that this has been done and q(x) is a pseudoscalar. The topological 
density has the classical continuum limit 

1 Ni d/2 

Illfo o^ 2 ^ = (47r) d / 2 (d/2)\ e v iUl -M2 u d/2 ^ F^iviF^.^ • • • Ffj, d/2 u d/2 (x) ■ (4.2) 

At the moment, we cannot prove the above conjecture in the non-perturbative 
level. However, we see that the conjecture holds to all orders in perturbation theory; 
the following theorem guarantees that a gauge invariant topological field is unique 
(up to a total divergence) under certain conditions: 

8 For G = U(l), the construction of ref. [29] can be generalized to arbitrary dimensions [26]. The 
equivalence of eq. (4.1) with eq. (2.11) for G = U(l) has been shown [26]. See also ref. [30]. 
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Theorem 4.1 Let p(x) be a local gauge invariant pseudoscalar field (which is trans- 
lational invariant) on the infinite lattice whose dependences on the lattice spacing a 
arise only though the gauge field. 9 If it is topological 

E S p( x ) = 0, (4-3) 

and the classical continuum limit lim a ^ p(x)/a d vanishes, then to all orders in per- 
turbation theory, 

p(x) = d;e,(x), (4.4) 

where i^x) is a local gauge invariant axial vector current. 

Proof. Our proof is rather similar to the cohomological argument of ref. [5]. We 
expand p(x) with respect to the bare gauge coupling constant go introduced by 
U(x,fjb) = e 9oA ^ x) 

oo 

p( x ) = 

k=l 

K - yi,—,Vk 

where A^(x) = A a ^(x)T a . 

First consider p^(x). Since p(x) is gauge invariant, p^(x) is invariant under 
the linearized gauge transformation 

A^x) ^ A^+d^x), (4.6) 

and also under the constant gauge transformation 

A^x) -> A^x) + [u, A^x)}. (4.7) 

Moreover, since p^\x) is a local topological pseudoscalar field and eq. (4.6) is the 
gauge transformation in abelian theory, one can invoke the cohomological analysis 
in abelian theory. The result is 

P (1 \x) = d*/;\x), £«(*) = g Q Y,$\xMK(v)- (4-8) 

y 

The local axial vector current tfi(x) is invariant under eqs. (4.6) and (4.7). A key 
observation is that, from £^\x), one can construct a field £^\x) such that it is 
invariant under the original non-abelian gauge transformation and its lowest-order 

9 Recall that in the classical continuum limit the gauge potential is introduced as U(x,n) = 
■Pexp[a J Q dt A^(x + (1 — t)afi)} where a is the lattice spacing. 
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O(go) term coincides with iff{x). This can be accomplished by substituting the 
gauge potential A a (y) in eq. (4.8) by the expression [5] 

Afa, y) = - tr{T a \l - W(x, y)U(y, »)W(x, y + A)" 1 ! }, (4.9) 

where W(x,y) is the ordered product of the link variables from y to x along the 
shortest path that goes first in direction 1, then direction 2, and so on. Note that 
A a {x, y) behaves gauge covariantly under the original non-abelian gauge transforma- 
tion. Thus the resulting expression, 

%\x) = 9 T,*i ) (x,y)l%(x,y), ( 4 -!0) 

y 

is invariant under the non-abelian gauge transformation due to the invariance of iff (x) 
under eq. (4.7). Moreover, since 

A„(x,y) = A^y) + dyu(x,y) + O(g ), (4.11) 

with u(x, y) the oriented line sum of the gauge potential from y to x, the invariance 
under eq. (4.6) implies that iff{x) = iff(x) + 0{gl). w Using iff{x), we may define 
a local gauge invariant pseudoscalar field 

p(x)-d*/ff(x). (4.12) 

which has identical properties with p(x) except that it starts with O(qq) term. Thus 
we can repeat the above argument from eq. (4.5) for the field (4.12). This time, 
however, the perturbation series analogous to eq. (4.5) starts from k = 2. 

In this way, we repeat the steps from eq. (4.5) to eq. (4.12) by eliminating the 
lowest-order term of the topological field until the first order term becomes 0(g^ 2 ); 
here a new situation arises. The cohomological analysis (with the fact that it is a 
pseudoscalar) tells that 

^w-^A^C^ + /*! + *)••• 

xKd/lvi/*(n + Ai + h H V /id/2-1 + Od/2-1) 

+d;tf 2 \ (4.i3) 

[F£ v (x) = d^Al(x) — d v A°^(x) denotes the linearized field strength] is a possible form 
of p(x) d*Jff (x) . However, since the continuum limit of p(x), lim a _^ p{x) / a 

vanishes, we infer that the constants c ai "' ad / 2 vanish, c ai "' ad ' 2 = 0. Thus we again 
have a total divergence. Further repeating the above procedure, we finally establish 

10 The current ftff(x) so constructed is not an axial vector under the lattice symmetries. However 
we can always enforce this by taking average over lattice symmetries. 
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Going back to eq. (4.1), we note that both A°°(x) and q(x) are a local gauge 
invariant topological pseudoscalar field (for the latter, those properties follow from 
the construction of q(x) [29]). Moreover, they have the same classical continuum 
limit (4.2). Thus, applying theorem 4.1 to A°°(x) — q(x), we see that the conjecture 
holds to all orders in perturbation theory. 

Now, in the proof of theorem 3.1 in abelian theory, every steps are valid even 
for non-abelian theories, except for the crucial relation (3.4), namely Yxer A°°(x) is 
an integer. With our conjecture 4.1 for non-abelian cases, this last condition is also 
satisfied; YxerQ( x ) is Liischer's topological charge on a periodic lattice which is an 
integer. So, repeating the proof for theorem 3.1, we have 

Corollary of conjecture 4.1 For general G, if the lattice is sufficiently large com- 
pared to the localization range g of the Dirac operator, say L/g>n, 



where fc M (x) is a gauge invariant periodic current on r. The current k^(x) moreover 
satisfies the bound 



with constants K\ and V\. The topological density q(x) is given by Liischer's topolog- 
ical density [29] and its higher dimensional extensions. 

This corollary states that the basic structure of axial anomaly on finite lattices is 
identical that on the infinite lattice. Summing eq. (4.14) over the lattice T, one has 
an equality between the index of the Dirac operator (2.3) and the geometrically- 
defined lattice topological charge [29]. This equivalence ("lattice index theorem") 
has been thought to be true for long time since the analyses in refs. [27, 31]. Our 
argument provides a further support for this equivalence. 

5. Conclusion 

In this paper, we have studied the axial anomaly defined on a finite-size lattice by 
using a Ginsparg- Wilson Dirac operator. For G = U(l), we show that the basic 
structure of axial anomaly on the infinite lattice, which has a quite analogous form 
to the continuum counterpart, persists even on a sufficiently large finite-size lattices. 
For general G, we conjectured that the axial anomaly on the infinite lattice is basically 
given by Liischer's topological density; actually this holds to all orders in perturbation 
theory. With this conjecture, we showed that this structure again persists even on 
finite-size lattices. Since Liischer's topological density is a geometrically natural 
definition of the Chern form in lattice gauge theory (note that it is proportional 
to strT ai • • -T ad / 2 ), our analysis indicates that the basic structure of axial anomaly 
in continuum theory is quite robust and it persists even in a system with finite 




(4.14) 




(4.15) 
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ultraviolet and infrared cutoffs. Of course, we indicated this persistency only in 
a framework with the Ginsparg- Wilson relation. To understand precise conditions 
on the formulation for this persistency to hold is an interesting open question; for 
example, one may enlarge the set of formulations by using the generalized Ginsparg- 
Wilson relation [32]. 

In the gauge invariant lattice formulation of abelian chiral gauge theories [15], 
a knowledge on the structure of U(l) gauge anomaly on finite-size lattices was of 
crucial importance. Recalling this fact, we believe that our analyses will be useful in 
extending the construction of ref. [15] to non-abelian gauge theories. 

H.S. would like to thank Takanori Fujiwara, Takahiro Fukui, Yoshio Kikukawa 
and Martin Luscher for valuable discussions. We are grateful to Kazuo Fujikawa for 
a careful reading of the manuscript. 
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